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Forewa rd 


The  Gulf  Stream  Meanders  Experiment  (GSME)  was  a  multi-phase" 
research  program  with  the  objective  of  obtaining  detailed  kinematical 
and  dynamical  descriptions  of  the  mesoscale  Gulf  Stream  meanders  which 
occur  along  the  continental  margin  of  the  southeastern  United  States. 

The  research  program  was  composed  of  a  theoretical  component,  which 
began  in  early  1977,  and  an  observational  component,  which  provided 
intensive  and  extensive  views  of  the  currents  in,  and  the  hydrographic 
structure  of  the  Gulf  Stream  between  Charleston,  South  Carolina,  and 
Cape  Hatteras,  North  Carolina,  during  1979. 

The  present  Technical  Report  describes  a  theoretical  study  of  sub- 
inertial  Gulf  Stream  fluctuations.  Thi3  was  one  contribution  to  the 
theoretical  component  of  the  GSME.  This  studv  was  undertaken  to  provide 
a  basic  description  of  the  structure  and  dispersion  properties  of  stable, 
Rossb' -like  waves  propagating  in  a  baroclinic,  laterally-sheared  western 
boundary  current.  With  a  continuously  stratified  mean  flow  and  variable 
bott<""r.  topography,  this  model  represents  a  significant  step  towards 
realism  in  tiie  theory  of  oceanic  long  wave  motions  along  a  continental 
margin.  The  results  of  this  study  indicate  the  importance  of  the  inter¬ 
play  among  the  density  field,  bottom  topography  and  mean  current  in  deter¬ 
mining  the  normal  modes  that  may  propagate  in  the  Gulf  Stream.  Taken 
together  with  field  observations  made  in  the  GSME,  the  theoretical  results 
suggest  that  the  strongly  sheared,  cyclonic  Gulf  Stream  frontal  zone  is  at 
the  heart  of  mesoscale  Gulf  Stream  variability.  Theoretical  studies  of 
spatially  growing  unstable  wave  motions  in  a  baroclinic  Gulf  Stream  are 
presently  underway  at  The  University  of  North  Carolina  to  broaden  our 
understanding  of  the  complex  dyna&nicai  nature  of  Gulf  Stream  meanders. 


Abstract 


The  effects  of  a  realistic  western  boundary  current  on  the  alongshore 
propagation  of  sub-inertial  waves  trapped  by  a  sloping  bottom  topography 
are  studied  using  a  numerical  model  incorporating  realistic  bottom  topo¬ 
graphy  and  a  current  field  which  is  in  thermal  wind  balance  with  the  den¬ 
sity  field.  This  models  the  Gulf  Stream  as  it  flows  along  the  continental 
slope  off  North  Carolina.  The  mean  state  velocity  and  density  fields  do 
not  vary  alongshore  and  are  continuous  in  the  horizontal  as  well  as  the 
vertical.  The  linearized,  inviscid  equations  of  motion  for  small  amplitude 
disturbances  yield  a  single  governing  equation  for  the  perturbation  pres¬ 
sure.  This  equation  is  solved  using  a  marching  method  for  elliptic  problems 
The  dispersion  relations  are  obtained  by  searching  for  the  resonance  respons 
of  the  system  to  an  arbitraty  uniform  forcing  term.  Four  discrete  stable 
modes  of  Rossby-like  waves  are  identified,  all  propagating  in  the  upstream 
direction.  A  mode-coupling  resonance  is  found  between  the  first  two  modes. 
For  small  wavenumber,  the  first  mode  is  trapped  within  the  frontal  zone  on 
the  cyclonic  side  of  the  mean  current  with  a  smaller  barotropic  component 
over  the  shelf.  This  "frontal  trapping"  is  due  to  the  quasi-geostrophic 
potential  vorticity  gradient  in  the  mean  current.  For  large  wavenumber, 
the  first  mode  becomes  primarily  a  barotropic  shelf  wave.  The  second  mode 
is  a  purely  barotropic  shelf  wave  for  small  wavenumbers,  but  becomes  a 
f rontal-trapped  wave  at  large  wavenumbers.  For  the  higher  wave  modes,  most 
of  the  energy  is  trapped  near  the  surface  in  the  frontal  zone. 
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1.  Introduction 


Much  research  in  physical  oceanography  during  the  past  few  years 
has  been  devoted  to  the  subject  of  coastal-trapped  waves.  This  gen¬ 
eral  classification  of  sub-inertial  frequency,  topographically  trapped 
waves  can  be  broken  down  into  several  subclassifications  under  certain 
limiting  conditions  (Wang  and  Mooers,  1976). 

In  the  simple  case  of  an  ocean  of  constant  depth  with  a  vertical 
wall  at  the  coast,  these  waves  are  coastal  Kelvin  waves.  In  a  homog¬ 
enous  ocean  over  a  constantly  sloping  bottom  topography,  the  general 
coastal -trapped  wave  reduces  to  the  quasi-geostrophic  edge  wave  (Reid, 
1958).  If  the  slope  of  the  bottom  topography  changes  away  from  the 
coast,  such  as  a  finite  width  continental  shelf  adjoining  a  constant 
depth  deep  sea  region,  the  familiar  barotropic  continental  shelf  wave 
results  (Robinson,  1964;  Mysak,  1967;  Buchwald  and  Adams,  1968). 

The  structure  of  the  coastal -trapped  wave  becomes  more  complicated 
in  the  presence  of  a  laterally  sheared  barotropic  mean  current  (Niiler 
and  Mysak,  1971;  Grimshaw,  1976;  Brooks  and  Mooers,  1977a, b)  or  in  the 
presence  of  a  level  stratification  over  a  sloping  bottom  topography 
(Allen,  1975;  Wang,  1975;  Wang  and  Mooers,  1976;  Clarke,  1977;  Huthnance 
1978).  Mysak  (1980)  presents  a  comprehensive  review  of  the  theory  of 
coastal-trapped  waves. 

This  study  extends  the  present  theory  to  include  the  effects  of 
a  continuously  stratified,  laterally  and  vertically  sheared  western 
boundary  current  flowing  along  the  continental  slope. 

Reid  (1958)  first  described  quasi-geostrophic  edge  waves  trapped 
over  a  constantly  sloping  continental  shelf  of  infinite  extent.  These 
low  frequency  waves  are  supported  by  the  potential  vorticity  gradient 
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associated  with  the  sloping  bottom.  They  are  right  bounded,  that  is, 
they  propagate  their  phase  with  the  shallow  water  to  their  right,  in 
the  northern  hemisphere.  Reid  called  these  waves  "quasi-geostrophic" 
because  the  balance  in  the  momentum  equations  was  essentially  geostrophic . 

Continental  shelf  waves  were  first  observed  by  Hamon  (1962)  as 
non-barometr ic  sea  level  disturbances  that  propagated  northward  along 
the  Australian  coast  from  Sidney  to  Coff's  Harbor.  Robinson  (1964) 
developed  a  simple  linear  theory  to  explain  Hamon 's  observations. 
Robinson's  model  incorporated  a  constantly  sloping  continental  shelf 
of  finite  width  adjoining  a  constant  depth  deep  sea  region.  He  called 
these  low-frequency,  vorticity  wave  motions  continental  shelf  waves 
because  they  were  trapped  over  the  shelf  by  the  sloping  bottom  and 
decayed  exponentially  away  from  the  shelf. 

Mysak  (1967,  1968)  showed  that  Reid's  quasi-geostrophic  edge  waves 
and  Robinson's  continental  shelf  waves  were  closely  related.  Both  waves 
are  special  cases  of  a  topographic  Rossby  wave.  At  short  wavelengths, 
the  shelf  waves  become  quasi-geostrophic  waves,  since  they  are  trapped 
more  closely  to  the  coast  and  cannot  "see"  the  edge  of  the  shelf. 

Mysak  (1967)  considered  the  effects  of  atmospheric  pressure  variations 
on  the  generation  of  shelf  waves.  He  also  included  a  simple  stratifica¬ 
tion  and  current  in  the  deep  sea  region  of  Robinson's  (1964)  model 
geometry. 

Buchwald  and  Adams  (1968)  extended  these  linear  theories  to  include 
a  more  realistic  bottom  topography.  Using  an  exponential  shelf  and 
slope  depth  profile,  they  showed  that  the  offshore  trapping  scale  for 
the  shelf  wave  modes  was  the  shelf  width.  Adams  and  Buchwald  (1969) 
in  the  same  model  geometry,  showed  that  wind  stress,  rather  than 


atmospheric  pressure,  was  responsible  for  the  generation  of  shelf  waves. 
Gill  and  Schumann  (1974)  gave  a  more  detailed  analysis  of  this  genera¬ 
ting  mechanism.  They  showed  that  the  alongshore  component  of  the  wind 
drove  an  likman  transport  in  the  surface  mixed  layer.  Since  shelf  waves 
arc  non-d i vergent ,  conservation  of  mass  requires  that  there  is  a  return 
flow  through  the  interior.  It  is  this  return  flow  over  the  strongly 
sloping  bottom  that  changes  the  vorticity  by  cortex  stretching,  and  henc 
generates  shelf  waves. 

Recent  theoretical  studies  have  shown  that  lateral  boundary  current 
shear  strongly  affects  the  propagation  of  barotropic  shelf  waves.  In 
areas  such  as  the  Gulf  Stream  off  the  east  coast  of  the  United  States, 
the  horizontal  shear  of  the  mean  current  can  be  comparable  to  the  Cor¬ 
iolis  parameter  f.  Shelf  waves  can  be  significantly  advected  by  the 
mean  current  and  have  their  propagation  characteristics  strongly  mod¬ 
ified.  1'hey  can  extract  energy  from  the  mean  current  through  the  pro¬ 
cess  of  barotropic  instability.  The  strong  shear  of  the  current  can 
support  a  new  class  of  shear  waves.  Niiler  and  Mysak  (1971)  examined 
the  trapped  wave  solutions  for  the  case  of  a  constantly  sloping  shelf 
of  finite  width  with  a  V-shaped  alongshore  mean  current.  This  topog¬ 
raphy  and  current  structure  was  chosen  to  model  the  Gulf  Stream  as  it 
flows  along  the  southeastern  coast  of  the  United  States.  They  found 
that  the  direction  of  shelf  wave  propagation  was  reversed  by  the  mean 
flow  at  a  high  wave  number  cut-off.  There  also  exists  a  class  of 
shear  waves  that  always  travel  in  the  downstream  direction.  The  dis¬ 
persion  curves  for  the  shelf  wave  and  the  shear  wave  cross  in  a  partic¬ 
ular  range  of  wavenumber,  where  the  waves  become  unstable.  At  wave- 
numbers  above  this  range,  the  two  dispersion  curves  change  families,  the 
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shelf  wave  curve  becoming  a  shear  wave,  and  the  shear  wave  curve  becom¬ 
ing  a  shelf  wave. 

Brooks  and  Mooers  (1977b)  studied  stable  barotropic  shelf  waves 
over  an  exponential  depth  profile  with  an  exponential  alongshore  mean 
current.  They  also  found  that  the  direction  of  shelf  wave  propagation 
was  reversed  by  the  mean  flow  at  a  high  wavenumber  cut-off.  Brooks 
(1975),  Grimshaw  (1976)  and  McKee  (1977)  discussed  the  problem  of  crit¬ 
ical  layers  in  the  mean  flow;  that  is,  a  level  in  the  mean  flow  where  the 
wave  phase  speed  equals  the  mean  current  velocity.  In  addition  to  the 
discrete  spectrum  of  shelf  wave  modes,  there  is  a  continuum  of  allow¬ 
able  critical  layer  solutions  in  the  dispersion  plane,  with  phase  speeds 
lying  in  the  range  of  the  current  speed. 

All  of  the  shelf  wave  models  discussed  above  have  dealt  with  a 
homogenous  ocean.  The  presence  of  stratification  can  alter  the  proper¬ 
ties  of  shelf  waves  and  also  can  allow  other  types  of  waves  to  exist. 

A  number  of  models  have  incorporated  a  level  stratification  in  the 
study  of  coastal -trapped  waves. 

Mysak  (196?)  and  Gill  and  Clarke  (1974)  examined  long  waves  in  a  two- 
layer  ocean  with  only  the  top  layer  extending  over  the  shelf.  It  was 
shown  by  Allen  (1975)  and  Wang  (1975)  that  the  barotropic  and  baroclinic 
modes  were  coupled  when  the  two- layer  stratification  extended  over  a 
sloping  continental  shelf.  Allen  (1975)  showed  that  the  strength  of 
this  coupling  depended  on  the  ratio  of  the  internal  radius  of  deforma¬ 
tion  to  the  length  scale  of  the  topographic  variation.  Both  Allen  (1975) 
and  Wang  (1975)  found  that  where  the  dispersion  curves  of  the  shelf 
waves  and  the  single  Kelvin  wave  appeared  to  cross,  there  was  actually 
a  change  of  modal  structure  of  the  waves.  Wang  (1975)  called  this 
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coupling  between  the  shelf  waves  and  the  internal  Kelvin  wave  a 
"resonance  coupling".  This  terminology  was  introduced  by  Eckart  (1962) 
in  a  study  of  internal  waves. 

Wang  ami  Moocrs  (1976)  also  found  this  resonance  phenomenon  in 
a  continuously  stratified  model.  They  found  that  in  the  limit  of  a 
vanishing  coastal  wall,  topographic  Rossby  waves  were  the  only  class  of 
sub- inertial  frequency,  coastal-trapped  wave  motion.  The  topographic 
Rossby  wave  reduces  to  a  barotropic  continental  shelf  wave  in  the  special 
case  of  small  stratification,  and  to  a  bottom-trapped  wave  (Rhinos, 

1970)  in  the  special  case  of  large  stratification.  Huthnance  (1978) 
showed  analytically  that  for  a  monotonical ly  increasing  depth  profile 
with  a  level  density  stratification  there  was  only  one  infinite  discrete 
set  of  trapped  sub- inertial  modes  with  frequency  decreasing  as  mode 
number  increased.  Clarke  (1977)  showed  that  in  the  long  wave  limit, 
these  waves  were  a  hybrid  between  shelf  waves  and  internal  Kelvin  waves. 
These  models ,  however,  have  not  considered  the  presence  of  sloping 
isopycnals  with  an  associated  vertically  sheared  mean  current  in  thermal 
wind  balance  with  the  mean  density  field. 

The  Gulf  Stream  region  off  the  southeastern  coast  of  the  United 
States  is  characterized  by  the  presence  of  strongly  sloping  density  sur¬ 
faces  combined  with  strong  lateral  and  vertical  current  shear.  Orlanski 
(1969),  using  a  two- layer  model,  and  Orlanski  and  Cox  (1973),  using  a 
continuous  stratification,  studied  barocl inical ly  unstable  waves  in  a 
western  boundary  current,  but  neglected  stable  waves. 

Ranc  and  llsueh  (1980)  investigated  the  role  of  a  density  front  in 
determining  the  dispersion  characteristics  of  stable  topographic  Rossby 
waves.  they  found  a  "complementary  mode"  in  which  a  barotropic  wave 
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over  the  shelf  matched  with  a  baroclinic  wave  at  the  density  front. 

Bane  (1980)  developed  a  two- layer  model  which  incorporates  a  general 
continental  shelf  and  slope  topography  with  a  thermal -wind  mean  current 
in  the  upper  layer.  He  identified  four  types  of  stable,  sub-inertial 
frequency,  llosshy-like  waves.  They  are  the  barotropic  shelf  wave,  the 
quas i-geostrophic  edge  wave,  the  complementary  mode  edge  wave,  and  a 
new  frontal -trapped  wave,  so-called  because  its  amplitude  is  trapped 
within  the  frontal  zone  on  the  cyclonic  side  of  the  mean  current.  The 
dispersion  curves  for  this  model  exhibit  the  mode-coupling  resonance 
phenomenon,  allowing  the  characteristics  of  the  component  waves  to  be 
interchanged  along  a  composite  dispersion  curve.  An  obvious  limiting 
feature  of  these  two-layer  models  is  that  they  lack  sufficient  vertical 
resolut ion. 

there  is  a  great  deal  of  evidence  to  support  the  existence  of 
shelf  waves  in  general  (liamon,  1962,  1966;  Mooers  and  Smith,  1968; 
Cartwright,  1969,  Cutchin  and  Smith,  1973)  and  particularly  off  the 
coast  of  North  Carolina  (Mysak  and  Hamon,  1969;  Brooks,  1978).  It  has 
been  suggested  that  mcandcrings  of  the  Gulf  Stream  in  this  area  may  be 
related  to  coastal -t rapped  waves  (Brooks,  1978;  Brooks  and  Bane,  1978; 
Bane  and  Brooks,  1979a).  Stable,  propagating,  wave-like  meander  patterns 
in  the  inshore  edge  of  the  Gulf  Stream  have  been  observed  to  be  a  dom¬ 
inant  mode  of  oscillation  of  the  Stream  (Legeckis,  1979).  The  existing 
theoretical  models,  however,  have  not  dealt  adequately  with  the  complex 
dynamics  found  in  this  region.  On  the  inshore,  or  cyclonic,  side  of 
the  Gulf  Stream,  both  the  strongly  sloping  density  surfaces  and  the 
large  horizontal  shear  of  the  mean  flow  contribute  in  the  same  sense  as 
the  sloping  bottom  to  the  ambient  potential  vorticity  gradient,  giving 
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rise  to  dual  trapping  mechanisms.  This  second  trapping  mechanism  is 
the  gradient  of  the  quasi-geostrophic  potential  vorticity  (cf,  Pedlosky, 
1970).  it  is  the  two  layer  analogue  of  the  quasi-geostrophic  potential 
vorticity  gradient  that  supports  the  frontal  trapped  waves  described 
by  Bane  (1980).  The  maximum  velocity  of  the  Stream  is  on  the  same 
order  as  the  phase  speed  of  topographic  Rossby  waves,  but  in  the  oppo¬ 
site  direction,  so  that  one  would  expect  significant  advection  of  the 
waves.  It  has  been  shown  that  when  both  stratification  and  topography 
arc  present,  the  barotropic  and  baroclinic  modes  are  coupled  (Allen, 
1975).  Thus,  the  presence  of  a  western  boundary  current  in  a  coastal- 
trapped  wave  model  can  alter  the  properties  of  the  coastal-trapped  wave 
modes,  and  can  allow  the  existence  of  a  new  class  of  frontal-trapped 
waves.  Moreover,  one  expects  the  coastal-trapped  waves  to  be  coupled 
with  the  frontal-trapped  waves  due  to  the  stratification.  To  aid  in 
the  study  of  these  effects,  the  present  model  includes  a  realistic 
western  boundary  current  flowing  along  a  continental  slope. 

Mysak  (1980)  points  out  that  the  presence  of  a  wide  flat  shelf, 
as  is  found  off  the  coast  of  the  Carolinas,  enhances  the  mode-coupling 
resonance  found  by  Wang  and  Mooers  (1976) .  Bane  (1980)  suggests  that 
this  is  due  to  changes  in  the  dispersion  properties  of  the  different 
component  waves  with  changing  shelf  width.  This  model  simulates  the 
Gulf  Stream  as  it  flows  along  the  continental  slope  off  Onslow  Bay, 

North  Carolina,  where  the  shelf  is  very  wide.  In  this  continuous  strat¬ 
ification  model,  the  different  wave  types  found  in  Bane's  (1980)  two- 
layer  model  are  more  subtly  coupled  and  are  more  difficult  to  distinguish 
however,  the  frontal-trapped  wave  is  clearly  present.  A  hybrid  between 
a  topographic  Rossby-likc  wave  and  a  frontal-trapped  wave  is  found  also. 


The  next  section  describes  the  formulation  of  the  governing  equa¬ 
tion  and  the  boundary  conditions.  Section  3  describes  the  mean  state 
for  the  model,  and  section  4  outlines  the  numerical  solution  scheme. 
The  results  of  the  computations  are  presented  in  section  5,  and  are 
discussed  in  section  6.  The  conclusions  are  summarized  in  section  7. 
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Formulation 


Consider  a  model  eastern  continental  shelf  and  slope  region  with 
a  steady  alongshore  flow  in  geostrophic  and  hydrostatic  balance  with 
the  density  field.  A  straight  coastline  is  assumed  along  the  y-axis, 
with  the  x-axis  pointing  offshore  and  the  z-axis  pointing  vertically 
upward.  The  depth,  h(x) ,  is  a  monotonically  increasing  function  of  x 
only.  A  schematic  of  this  idealized  geometry  is  shown  in  Fig.  2.1. 

The  fluid  is  assumed  to  be  inviscid,  incompressible,  non-diffusive,  and 
continuously  density-stratified.  Variations  in  the  Coriolis  parameter, 
f,  are  neglected.  The  equations  of  motion  for  small  amplitude  disturb¬ 
ances  are 
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where  subscripts  denote  differentiation,  and 
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Following  the  notation  of  Mooers  (1975a, b),  the  dependent  vari¬ 
ables  are  separated  into  mean  and  fluctuating  components  (denoted  by 
overbars  and  primes,  respectively) : 
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v  =  v  +  v '  ,  u  =  u',  w  =  w', 

p  =  P0  +  p  +  p '  , 

p  =  p0  +  p  +  p'  . 

where  pQ  =  -pQgz 
and  (  P,  p')«  P0  • 

The  buoyancy  and  the  nonhydrostatic  pressure  are  defined  as 

b  =  -  ~  =  b  +  b  +  b' 

P0  o 

1*  - 

it  =  -  =  it  +  n ' 

Po 

The  equations  of  motion  are  time  averaged  and  the  wave  motion  is  sep¬ 
arated  from  the  mean  flow  (see  Appendix  Al) .  Reynolds  stresses  and 
non-linear  or  wave-wave  interaction  terms  are  neglected  and  a  Boussinesq 
approximation  is  adopted. 

The  mean  flow  is  governed  by 

fv  =  Sx  (2.1) 

0  =  -  tiz  *  b  (2.2) 

Together,  (2.1)  and  (2.2)  yield  the  thermal  wind  relation 

f  v,  -- 

2 

where  M  =  b^  is  the  horizontal  analogue  to  the  Brunt-Vaisala  frequency, 
N2  =  *&>z  =  b  . 

Po  2 


A  hydrostatic  balance  is  assumed,  since  only  sub- inertial  frequency 
motions  are  being  considered.  The  wave  motion,  without  the  primes, 
is  governed  by 


+  vu  -  fv  =  -IT 
t  y  x 


v  +  vv  +  (f+v  )u  +  v  w  =  -7T 
t  y  x'  z  y 


0  =  -  n  ^  +  b 


(2.5) 

(2.4) 

(2.5) 


u  +  v  +  w  =  0 
x  y  z 


(2.6) 


?  ' 

b  +  vb  +  M“u  +  N“w  =  0 

t  y 


(2.7) 


2 

Because  ot  the  strong  horizontal  density  gradients,  M  must  be  included 

in  (2.7);  therefore,  from  the  thermal  wind  relation,  v^  must  be  retained 

in  (2.4).  The  strong  lateral  shear  of  the  Gulf  Stream,  which  augments 

the  background  potent ia 1  vorticity  of  the  region,  requires  the  inclu- 

s ion  of  v  in  (2.4)  . 
x 

Only  alongshore  propagating  disturbances  are  considered;  therefore, 
all  perturbation  variables  (u,v,w,n,b)  are  assumed  to  take  the  form 

ci(o0t  +  iy)  their  (y,t)  dependence,  where  o0  is  the  wave  fre- 

qucncy,  t  is  the  alongshore  wavenumber,  and  i  =  (-1)  '  .  Differentia¬ 
tion  with  respect  to  t  and  y  thus  becomes  equivalent  to  multiplcation 
by  io  and  If,  respectively.  The  intrinsic  frequency  is  the  frequency 
seen  by  a  fluid  particle  at  a  particular  x  and  z,  and  is  defined  as 
<i(x,z)  =  o  +  v(x,z)£  .  The  mean  current,  v,  is  always  positive  (a 
northward  flowing  current),  so  that  for  a  wave  propagating  in  the 
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positive  y  (downstream)  direction  (either  aQ  or  l  negative,  but  not 
both),  there  exists  the  possibility  of  a  critical  surface  where  a(x,z)  = 
0,  i  .e. ,  the  wave  phase  is  propagating  at  the  same  speed  as  the  local 
medium.  For  this  study,  only  stable  waves  propagating  in  the  negative 
y  (upstream)  direction  will  be  considered,  so  that  c>o  and  l  are  always 
positive,  and  critical  levels  are  avoided.  Equations  (2.3)  through 
(2.7)  now  become 


iou 

-  fv  =  -It 

(2.8) 

X 

iuv 

+  (f+v  )u  +  v  w  =  -i£n 

(2.9) 

X  z 

0  = 

-tt_  +  b 

(2.10) 

u  + 

X 

ifv  +  w  =  0 

(2.11) 

7  l 

iob  + 

M“u  +  N“w  -  0 

(2.12) 

Eliminating  the  buoyancy,  b,  between  (2.10)  and  (2.12),  then 
solving  (2.8)  through  (2.10)  for  u,  v  and  w  in  terms  of  n,  gives 


u  - 

K 

i 

z 

r-v 

=4 

X 

..2 

+  oM 

z  ' 

(2. 

.13) 

V  = 

M 

7  2  2 
\  N  °f  - 

M1) 

+  of^N^T^ 

(2. 

.14) 

w  = 

i 

A 

-  a  (a 

^  -  o~  +  ^fM^nj 

(2. 

.15) 

2/  2  2  \  4  2 

where  A  =  N  -  a  j  -  M  ,  and  =  f(f  +  v^)  is  the  square  of  the 
effective  local  inertial  frequency  (Mooers ,  1975a).  Substituting  (2.13) 
through  (2.15)  into  the  continuity  equation  (2.11)  and  applying  the 
thermal  wind  relation  yields  the  governing  equation  for  the  perturbation 


pressure,  n : 
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if  f MJ  DA 
aA  Dz 


=  0 


(2.16) 


The  boundary  conditions  chosen  for  sub-intertial  motion  are  as 
follows.  At  the  surface  (z  =  0),  a  rigid  lid  is  imposed,  i.e.,  w  =  0. 

The  rigiil  lid  approximation  suppresses  any  wave-associated  sea  level 
elevations  and  free  surface  divergences  while  retaining  the  pressure 
gradients  along  the  surface.  This  has  the  effect  of  filtering  out  any 
external  gravity  waves.  Huthnance  (1978)  showed  that  the  free  surface 
divergence  effects  were  generally  very  small  for  sub-intertial  waves, 
and  that  the  eigenfrequencies  of  these  waves  were  increased  only  slightly 
by  invoking  the  rigid  1 i <1  approximation.  Using  (2.15),  the  boundary 
condition  at  z  =  0  hecomcs 


N,2\  '  (°f  ' 

At  z  =  -hfx),  there  is  no  flow  normal  to  the  bottom,  which  implies 
that  w  =  -uh^ ,  or  from  (2.13)  and  (2.15), 

2  2  if  2  ">  2 

(M  -  N  h  )  (tt  +  —  tt )  ♦  (M  h  -  (o;  -  o  ))  it 
X  X  o  x  f  z 


2<\ 

a  j  n2 


£f 

0 


M2h 


=  0  at  z=0. 


(2.17) 


=  0  at  z--h(x) .  (2.18) 


14 


At  the  coastline,  there  is  no  flow  normal  to  shore;  that  is,  at  x  =  0, 
u  =  0,  or  from  (2.13) 

M2  if  n  r  n  n  icn 

it  -  -r  u  +  —  ti  =  0  at  x=0.  (2.19) 

x  ..Z  Z  0 


Trapped  wave  solutions  arc  sought,  requiring  that  all  solutions  decay 
far  offshore;  that  is, 

tt  0  as  x  *  <*■  .  (2.20) 

Liquations  (2.16)  through  (2.20)  form  an  eigenvalue  problem  in 
and  l.  In  the  limiting  case  of  a  level  stratification  and  no  mean 
current  (v  =  0  and  =  0),  (2.16)  through  (2.20)  reduce  to  the  equa¬ 

tions  used  by  Kang  and  Moocrs  (1976)  and  by  Huthnance  (1978).  Equation 

2  2 

(2.16)  is  elliptic  for  the  values  of  N  ,  M  and  v^  encountered  in  a 
typical  western  boundary  current  profile  provided  c»0  <  f.  It  is  con¬ 
venient  to  solve  (2.16)  -  (2.20)  in  a  rectangular  domain;  therefore, 
the  independent  variables  (x,z)  are  transformed  to  (x,c)  ,  where  t;  = 
zll/h(x)  ,  II  being  the  maximum  depth  and  h(x)  the  depth  at  offshore 
distance  x.  The  governing  equation  (2.16)  then  becomes 


(2.21) 
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The  boundary  conditions  become 


at  r,  =  -H.  (2.24) 


The  variables  arc  made  dimensionless  by  the  relations 
*  x/L,  z*  =  z/ll ,  r.*  =  c/H,  h*  =  h/ii,  o*  =  o/f,  =  h,  a*  =  of/f. 


?*  ■>  7  7 

N“  -  N“/N-  ,  and  M“  =  M  /M  . 

max  max 


In  dimensionless  form,  equation  (2.21)  becomes, with  the  asterisks  dropped, 
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3.  Description  of  the  Model  Geometry 

The  background  mean  current  in  this  study  reasonably  approximates 

the  Gulf  Stream  as  it  flows  along  the  coast  of  North  Carolina 

(Richardson,  Schmitzand  Niiler,  1969).  It  is  essentially  the  same 

as  that  used  by  Orlanski  and  Cox  (1973)  as  their  initial  conditions. 

2  2- 

This  form  is  chosen  because  the  values  of  N  ,  M  ,  v  and  their  deriva¬ 
tives  arc  easily  evaluated  at  any  point  in  the  solution  domain. 

The  mean  velocity,  v,  is  assumed  to  be  in  geostrophic  balance  with 
the  density  field  which  is  determined  entirely  by  the  temperature.  The 
mean  velocity  is  in  the  positive  y  (northward)  direction,  and  docs  not 
vary  with  y  or  t.  The  analytic  expression  for  the  temperature  field 
used  by  Orlanski  and  Cox  (1973)  is 

-  x  2  f  z 

T ( x , ; )  =  T  -  AT(3-2(6x  +  l)c  '(rz-l)e  (3.1) 

where  the  parameters  T  ,  AT,  6  and  c  are  adjusted  to  give  the  desired 
temperature  profile.  The  density  is  given  by 

Ofx.z)  =  »■  ( i -oiT  (  x  ,  z) )  (3.2) 

where  «  is  the  thermal  expansivity  of  sea  water,  and  is  taken  to  be 
a  constant.  By  using  (5.1),  (3.2)  and  the  thermal  wind  relation,  the 
expression  for  the  mean  velocity  is  found  to  be 

get  AT  .2  - 6x  .  „  ,  ,  .  2ez 

-  fe. —  6  xc  (a/2  -  ez)c 

tr 


v  ( \  ,  z ) 


(3.3) 


IK 


The  horizontal  shear  of  the  mean  velocity  is 


-  .  ,  guAi  .2,,  -6x  ,,  2cz 

V  (X,z)  =  --T—  (5  (l-tfx)c  (.">/2-EZ)e 

X  T  t 


From  the  definitions  of  Brunt- Vaisala  frequency,  NT,  and  its  horizontal 
analogue,  M“ ,  it  is  found  that 


\“(x,z)  =  -guAT(3-2(<S.\+l)e  )  (2ez-l)ee‘ 


2  ">  -  fix 

M  (x,z)  =  -  2gcxAT(  6~xc  )(ez-l)e 


(3.6) 


The  bottom  topography  used  in  the  model  has  a  hyperbolic  tangent 


II  -11 

h(x)  =  !lwal,  +  — 0  [1  +  tanh  (s  (x-a) )  ] 


where  II  , ,  is  the  depth  at  the  coastline,  H  is  the  offshore  maximum 
wa  11  1  o 

depth,  a  is  the  offshore  distance  of  the  maximum  slope  and  s  is  propor¬ 
tional  to  the  maximum  slope. 

The  parameters  T  ,  AT,  <$ ,  and  e  for  the  current  and  density  fields, 
and  a,  s,  and  II  for  the  bottom  topography,  are  chosen  to  model  the 

t'.ulf  Stream  as  it  flows  along  the  continental  slope  off  Onslow  Bay,  Nor  h 
Carolina,  at  about  33.5°N  latitude  (Fig.  3.1).  A  current  section  from 
the  area,  after  Richardson  ct  al .  (1969),  is  shown  in  Figure  3.2.  For 
the  model  ,  the  deep  water  temperature  Tq  is  chosen  to  be  4°C.  The  ver¬ 
tical  temperature  difference  AT,  is  7 °C.  To  give  a  current  of  about 
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200  km  wide,  <5,  the  inverse  of  the  horizontal  scale  length,  is  set  equal 
to  3  x  ID' S  m  *.  The  vertical  scale  of  the  current  is  determined  by 
i ,  which  for  a  1500m  deep  current  is  equal  to  1,3  x  10  ^  u  The  in¬ 
shore  edge  of  the  current,  or  the  front,  is  placed  100  km  from  shore, 
which  is  just  over  the  shelf  break.  The  position  of  maximum  bottom 
slope,  a.  is  1-15  km  from  shore  and  the  maximum  slope  is  equal  to  0.05 
which  corresponds  to  s  -  2.5  *  10  The  resulting  velocity  and  temp¬ 
erature  fields  are  shewn  in  Figures  3.3  and  3.4  respectively.  The  depth 
of  the  coastal  wall  is  100  m,  or  2.5"  of  the  maximum  depth  H  =  4000m. 

The  offshore  coordinate  is  made  dimensionless  by  the  width  of  the  cur¬ 
rent,  1.  -  200  km,  and  the  vertical  coordinate  by  the  maximum  depth  Hq. 

The  offshore  boundary  is  imposed  at  a  distance  x  =  2L.  It  will  be  seen 
later  that  this  is  sufficiently  far  offshore,  as  all  solutions  die  out 
far  inshore  of  this  boundary.  The  presence  of  the  coastal  wall  is 
somewhat  unrealistic  and  is  purely  for  mathematical  convenience.  If 
the  depth  at  the  coast  is  allowed  to  go  to  zero,  some  of  the  coefficients 
in  (2.25)  become  infinite.  The  implications  of  a  finite  coastal  wall  will 
be  discussed  in  a  later  section. 


-1 .  Method  of  Solution 


The  governing  equation  (2.2S)  is  a  two  dimensional,  second  order, 
partial  differential  equation  in  a  rectangular  domain.  Solutions  to 
this  equation  subject  to  the  boundary  conditions  (2.26j  -  (2.29)  are 
obtained  for  a  particular  and  £  using  the  marching  method  for  elliptic 
partial  differential  equations  of  Roache  (1978).  For  a  fixed  l,  the 
eigenvalues  of  o  arc  found  by  searching  for  the  resonance  response  of 
the  system  to  a  progressive  wave  forcing  term  of  unit  amplitude. 

The  marching  method  is  "direct";  that  is,  it  produces  an  alge¬ 
braically  exact  answer  in  a  finite  number  of  steps,  as  opposed  to  an 
iterative  method  in  which  each  iteration  gives  a  better  approximation 
to  the  exact  solution.  To  solve  the  boundary  value  problem  for  a  fixed 
and  £,  the  solution  is  guessed  at  one  boundary  and  the  finite  dif¬ 
ference  form  of  the  equation  is  "marched"  through  to  the  opposite 
boundary.  The  results  of  the  march  are  then  compared  to  the  boundary 
condition  and  the  initial  guess  is  corrected  accordingly.  A  second 
march  then  gives  the  correct  solution.  A  complete  description  of  the 
algorithm  used  can  be  found  in  Roache  (1978). 

The  stability  of  the  marching  method  is  highly  dependent  on  the 
number  of  grid  points  in  the  direction  of  march,  J.  The  method  becomes 
unstable  for  large  .1,  primarily  due  to  machine  round-off  error  (see 
Roache,  1978) .  This  behavior  places  an  upper  bound  on  the  resolution 
of  the  method.  For  this  problem,  suitable  results  are  obtained  using 
lb  grid  points  in  the  £  direction,  the  direction  of  march,  and  31  points 
in  the  x  direction  (Figure  4.1).  The  use  of  a  larger  value  for  J  results 
in  unacceptable  error  propagation  through  the  march,  even  when  using 
double  precision  on  the  CRAY1  computer.  The  grid  spacing  used  gives 
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dimensionless  horizontal  and  vertical  resolutions  of  0.0667.  In  dimen¬ 
sional  coordinates,  this  gives  a  horizontal  resolution  of  13  1/3  km, 
since  the  offshore  boundary  is  placed  at  x  =  2L  =  400  km.  The  vertical 
resolution  is  a  function  of  x,  and  in  dimensional  coordinates,  Az  = 

0.0667  h ( x )  . 

The  dispersion  relations  are  found  by  repeatedly  solving  the  boun¬ 
dary  value  problem  for  different  values  of  oq  and  i.  Dimensionless 
wavenumber  f  is  set  at  40  increments  between  0  and  5.0.  For  each  value 
of  £,  the  boundary  value  problem  is  solved  for  60  values  of  oq  between 
0  and  0.6.  The  total  kinetic  energy  over  the  entire  solution  grid  is 
calculated  for  each  value  of  and  stored.  The  array  of  kinetic  energy 
values  for  that  particular  l  is  then  searched  to  locate  the  local  maxima. 
These  maxima  arc  considered  an  indication  of  a  resonance  response  of 
the  system  and  therefore  the  points  (aQ,  l)  corresponding  to  these  max¬ 
ima  lie  on  or  very  near  a  dispersion  curve.  Fig.  4.2  shows  a  plot  of 
integrated  kinetic  energy  versus  oq  for  l  =  2.0.  The  dispersion  rela¬ 
tion  is  obtained  by  contouring  the  resonance  points.  To  get  sufficient 
resolution  in  f requency-wavenumber  space  requires  the  solution  of  the 
boundary  value  problem  2400  times.  Each  solution  requires  about  one 
second  on  the  CRAYl ,  so  that  to  produce  an  entire  dispersion  relation 
requires  forty  minutes  of  CRAY1  CPU  time.  For  wavenumbers  greater  than 
5.0  the  horizontal  scales  of  the  motions  approach  the  horizontal  resolu¬ 
tion  of  the  numerical  method;  therefore  the  dispersion  relation  could 
not  be  extended  to  higher  wavenumbers. 
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This  method  is  tested  using  the  model  of  Wang  and  Mooers  (1976) 
for  a  flat  stratification.  Using  their  depth  profile 


h(x) 


O.00S  +  0.36  x  , 


-0.33  +  1.33  x 

1  .0 


0.42  -  x  -  0 

1.02  >  x  >  0.4 
x  >  1.0 


and  static  stability 

N2(z)  =  -(z  -  0 . 05) / (2  -  0.3) 2 


and  a  stratification  parameter  S  of  order  unity,  this  method  reproduces 
their  results  (Figure  4.3).  Assuming  that  the  results  of  Wang  and 
Mooers  (1976)  are  correct,  it  seems  safe  to  assume  that  this  method 
produces  sound  results. 


5.  Results 


The  dispersion  curves  obtained  for  the  model  are  shown  in  Figure 
5.1.  Five  curves  arc  distinctly  identifiable  from  the  eigenvalue 
search.  A  mode  numbering  convention  is  used  to  identify  the  curves  in 
the  dispersion  diagram.  The  uppermost  curve  is  called  "mode  0",  the 
next  highest  curve  is  called  "mode  1"  and  so  on.  The  mode  numbers  do 
not  imply  any  particular  wave  structure.  The  uppermost  curve,  mode  0, 
is  an  internal  Kelvin  wave  trapped  near  the  bottom  against  the  vertical 
coastal  wall.  All  of  the  energy  is  contained  in  the  alongshore  and 
vertical  velocities  This  mode  is  non-dispersive  and  for  wavenumber  l 
greater  than  1.0,  it  is  no  longer  in  the  range  of  sub-inertial  motion. 

The  other  modes  are  vorticity  controlled,  Rossby-like  waves  and 
are  supported  hy  two  wave  guides  -  the  sloping  bottom  topography  and 
the  sloping  isopyonals  and  associated  shear  of  the  mean  current.  The 
pressure  (n)  structure  for  mode  1  for  wavenumber  l  =  1.25  is  shown  in 
Figure  5.2(a).  Figures  5.2(b)  and  (c)  show  the  horizontal  velocities 
associated  with  this  mode.  In  the  wavenumber  range  0  <  t  <  1.5,  this 
wave,  consisting  of  a  barotropic  portion  trapped  over  the  shelf  which 
co-osc i  1  lutes  with  a  baroclinic  portion  trapped  within  the  frontal 
zone,  i.e.,  within  the  region  of  strongly  sloping  isopycnals  on  the 
cyclonic  side  of  the  mean  current,  is  a  hybrid  between  a  shelf  wave 
and  an  internal  Rossby-like  wave.  The  hybrid  nature  of  the  wave  de¬ 
creases  with  increasing  wavenumber,  that  is,  the  barotropic  component 
over  the  shelf  diminishes.  For  intermediate  wavenumbers  (1.5<f<3.5), 
this  mode  becomes  a  "fronta 1 -t rapped"  wave  with  most  of  the  energy 
trapped  within  the  frontal  zone  (Figure  5.3).  As  the  wavenumber  in¬ 
creases  still  further,  the  hybrid  structure  over  the  shelf  returns  and 
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the  frontal -trapped  wave  begins  to  die  out.  For  wavenumbei  <>3.75, 
the  frontal-trapped  structure  is  completely  gone,  and  this  mode  becomes 
a  topographic  Rossby  wave  (Figure  5.4). 

Mode  2  for  low  wavenumbers  is  a  topographic  Rossby  wave  over  the 
shelf  that  dies  out  in  the  frontal  zone.  The  pressure  and  velocity  struc¬ 
tures  for  this  mode  are  shown  in  Figure  5.5.  For  higher  wavenumbers 
(•£>3.75)  this  mode  begins  to  exhibit  baroclinic  structure  in  the  frontal 
zone  and  becomes  a  hybrid  wave  (Figure  5.6).  It  therefore  appears  that  modes 
1  and  2  exchange  characteristics  around  l  =  3.75.  This  is  similar  to  the 
resonance  coupling  found  by  Bane  (1980)  in  a  two-layer  Gulf  Stream  model. 

Mode  5  is  a  hybrid  wave  for  all  wavenumbers  considered  (Figure  5.7), 
while  mode  4  exhibits  the  frontal-trapped  structure  (Figure  5.8). 

Aside  from  the  single  internal  Kelvin  wave  mode,  there  appear  to 
he  three  distinctly  different  wave  structures.  The  first  is  a  topo¬ 
graphic  Rossby  wave,  which  is  confined  to  the  shelf  region  (i.e.,  a 
continental  shelf  wave)  and  dies  out  in  the  frontal  zone  of  the  mean 
current.  The  second  is  the  f ronta 1  -  trapped  wave,  which  is  essentially 
an  internal  Rossby  wave  supported  by  the  vorticity  gradient  resulting 
from  the  sloping  density  surfaces  and  shear  associated  with  the  mean 
current.  Over  the  shelf,  where  the  stratification  is  level,  there  is 
very  little  motion  associated  with  this  wave.  The  third  type  of  wave 
is  a  hybrid  between  the  first  two  and  has  components  over  the  shelf 
and  in  the  current  region  that  arc  of  comparable  magnitude.  These 
three  structures  can  be  understood  by  considerine  the  potential  vorticity 
control  mechanisms  which  support  the  wave  motion.  For  barotropic  motions 


25 

over  a  sloping  bottom,  conservation  of  potential  vorticity  dictates 
that 


(5.1) 


where  £  =  vx  _  u  is  the  vertical  component  of  the  relative  vorticity 
of  the  wave  motion.  This  is  the  case  over  the  shelf;  therefore,  the 
wave  vorticity  can  be  balanced  only  by  the  topographic  background  vor¬ 
ticity  f/h  and  only  topographic  Rossby  waves  are  allowed.  The  effect 
of  a  level  stratification  is  to  concentrate  this  background  vorticity 
gradient  in  the  lower  layers,  leading  to  bottom  trapping  of  the  wave 
motion,  as  in  Figure  5.5. 

The  presence  of  a  mean  current  with  both  lateral  and  vertical 
shear  gives  an  additional  vorticity  constraint,  the  conservation  of 
quasi -geostrophic  potential  vorticity  (cjf,  Pedlosky,  1979) 

m  =  iff  [  <«*«>  '  s;  ]  = 0  l=-2> 

A  _  A 

where  w  =  V *(ui  +  (v+v) j+wk)  is  the  total  vorticity  vector.  The  wave 
vurticity  can  now  bo  balanced  by  the  background  quasi-geostrophic 
potential  vorticity 

Vo 

Q  =  (fk  -  v  i)  . 


(5.3) 
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(whore  f  =  (f+v^)),  or  from  the  thermal  wind  relation 


1  /  4  2  2  \ 

Q  =  k  (  M  '  N  °f  ) 


(5.4) 


This  quantity  is  plotted  in  Figure  5.9.  The  quasi-geostrophic  poten¬ 
tial  vorticity  balance  is  then 


r«(w)* 


0 


p  p 

z  x 

where  Q'  =  f,  —  +  (w  -v  ) —  is  the  wave  quasi-geostrophic  vorticity. 

p0  y  2  po 

The  current  alone  can  support  internal  Rossby  waves  due  to  the  gradient 
in  Q.  This  gradient  is  strongest  near  the  surface  in  the  frontal  zone; 
hence,  the  wave  motion  may  be  concentrated  there.  The  strong  density 
stratification  in  the  current  effectively  isolates  these  waves  from 
the  bottom  topography  for  all  but  the  very  longest  wavelengths.  This 
can  be  seen  by  comparing  Figures  5.2(a)  and  Figure  5.3.  Mode  1  shows 
a  strong  hot  tom-trapped  component  at  l  =  1.25  (Figure  5.2(a))  which 
disappears  at  f  -  2.0  (Figure  5.3). 

For  a  level,  two  layer  stratification  over  bottom  topography, 

Allen  (1975)  has  shown  that  t.ie  strength  of  the  coupling  between  baro- 
tropic  and  baroclinic  modes  depends  on  a  parameter  A  =  ri/ 6  , 

D 

NH 

where  r.  =  j-  is  the  internal  Rossby  radius  of  deformation  and  =  H/H^ 
is  a  length  scale  for  the  bottom  topography.  Over  the  shelf  in  the 
present  model ,  N  ~  4  x  10  ^  sec  ^ ,  II  10^m,  f  “  10  ^sec  ^  and 
A  ~  100  km,  so  that  r.  ~  4km  and  A  ~  0.004.  In  the  current,  however, 

n  1 

N  ~  5  x  io'S  sec  1 ,  II  1000  m  and  r^  ~  50  km,  giving  A  =  0.5.  This 


-1 
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parameter  was  derived  for  a  level  stratification  and  is  not  totally 
applicable  here.  As  shown  above,  the  presence  of  sloping  isopycnals 
permits  the  frontal-trapped  wave  motion  that  does  not  exist  in  a  level 
stratification.  Coupling  between  this  type  of  motion  and  waves  supported 
by  the  sloping  topography  is  another  matter.  If  the  phase  speeds  of 
the  frontal-trapped  wave  and  the  topographic  Rossby  wave  are  similar 
for  a  given  wavelength,  the  two  waves  can  co-osci Hate,  resulting  in 
a  hybrid  wave.  Otherwise,  only  one  type  of  wave  may  exist. 

The  phase  speeds  for  modes  1,  2  and  3  at  low  wavenumbers  are  very 
close  to  the  first  three  barotropic  shelf  wave  speeds  for  the  same  depth 
profile.  This  is  to  be  expected,  since  for  very  long  waves,  — ,  u' 
and  w'  in  (2.3)  through  (2.7)  become  small.  Each  mode  (excluding  the 
internal  Kelvin  wave  mode)  has  a  high  frequency  cutoff,  resulting  in 
a  zero  group  speed  point  for  each  mode.  Table  5.1  gives  the  values 
of  the  cutoff  frequency  oc ,  the  corresponding  wavenumber,  and  the  phase 
speed  associated  with  the  zero  group  speed  of  each  mode  for  the  first 
three  Rossby- like  wave  modes.  The  higher-order  wave  modes  are  not 
adequately  resolved  by  the  solution  method. 


Table  5.1.  Values  for  high  frequency  cutoff,  for  the  first  three 
Rossby -like  wave  modes. 


Mode 

0 

c 

Dimens ional 
Frequency 

t 

Dimensional 

Wavelength 

Dimensional 
Phase  Speed 

1 

0.321 

2.62  x  10  ‘’sec  * 

2.25 

558.5  km 

200  km/day 

0.223 

1 .804  x  10  5scc  * 

3.7 

339.6  km 

84  km/day 

3 

0.15 

1.21  x  io  ^scc  * 

4.13 

304 . 3  km 

50  km/day 

0.  Discussion 


One  usually  thinks  of  a  wave  mode  in  terms  of  a  specific  modal 
structure;  for  instance,  the  number  of  nodes  in  the  horizontal  or 
vertical  pressure  structure.  This  concept  of  modal  structure  does 
not  seem  to  apply  to  the  modes  found  in  the  present  model.  Each  mode 
is  a  composite  of  several  different  wave  structures,  with  a  particular 
mode  exhibiting  different  structures  in  different  regions  of  the  dis¬ 
persion  plane.  Modes  1  and  2  show  a  coupling  resonance  phenomenon, 
wherein  the  two  modes  exchange  characteristics.  The  term  "resonance", 
which  was  first  used  in  this  context  by  Eckart  (1962)  ,  is  somewhat 
misleading  since  there  is  no  infinite  growth  in  the  wave  amplitude 
(cf  Allen,  1975)  . 

Similar  phenomena  have  been  observed  in  other  coastal-trapped  wave 
models.  Niilcr  and  Mysak  (1971)  found  that  the  dispersion  curves 
for  their  shelf  wave  modes  and  shear  wave  modes  change  families  where 
the  curves  appear  to  cross.  Allen  (1975),  Wang  (1975)  and  Wang  and 
Mooers  (1976)  found  similar  coupling  between  topographic  Rossby  wave 
modes  and  internal  Kelvin  wave  modes  in  models  that  included  a  level 
stratification  over  a  sloping  topography.  Bane  (1980),  in  a  two-layer 
Gulf  Stream  model,  found  that  the  dispersion  curves  for  a  particular 
topography/density/current  setting  may  be  interpreted  as  a  composite 
of  the  families  of  dispersion  curves  of  the  four  different  wave  types 
present  in  the  model.  These  four  wave  types  are  the  barotropic  shelf 
waves,  the  barotropic  quasigeostrophic  edge  waves,  the  complementary 
mode  edge  waves  and  the  frontal -trapped  waves.  In  the  present  contin¬ 
uously  stratified  model,  these  four  types  of  waves  are  not  as  easily 
distinguished  as  they  are  in  Bane's  (1980)  model.  This  is  due  to  the 
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fundamental  difference  between  a  two-layer  stratification  and  a  con¬ 
tinuous  stratification.  In  the  two-layer  case,  all  the  density  infor¬ 
mation  is  contained  in  a  single  discontinuity;  whereas  in  the  continu¬ 
ously  stratified  case,  the  density  information  is  distributed  over 
the  entire  water  column,  giving  much  greater  vertical  resolution.  Thus, 
the  interactions  between  stratification  and  topography  are  much  more 
subtle.  It  is  important  to  note  that  one  cannot  speak  of  the  mean 
velocity  and  mean  density  fields  independently,  as  they  are  intimately 
coupled  through  the  thermal  wind  relation.  In  this  particular  model, 
the  mean  density  field  was  specified  and  the  mean  velocity  field  was 
computed  from  it;  however,  one  could  have  just  as  easily  specified  the 
mean  velocity  field  and  then  computed  the  associated  density  field. 
Therefore,  when  we  speak  of  the  mean  current  or  the  mean  stratifica¬ 
tion,  we  are  referring  to  both  the  mean  velocity  and  density  fields. 

The  barotropic  shelf  wave  dispersion  curves  were  computed  for  the 
model  bottom  topography  and  arc  shown  in  Figure  6.1.  Comparing  these 
curves  with  those  of  Figure  5.1,  it  is  clear  that  the  presence  of  the 
mean  current  significantly  alters  the  propagation  of  the  shelf  waves. 

For  very  long  wavelengths  (i  ^  0.5)  both  sets  of  curves  are  roughly 
coincident.  Hi  is  is  to  be  expected,  since  for  very  long  wavelengths 


-  ,u,  and  w 

3y 

other  terms, 
curves  for  the 
curves.  It 
t  rans i t ion  to 


in  (2.5)  -  (2.7)  become  very  small  in  relation  to  the 
lor  smaller  wavelengths  (larger  wavenumber)  the  dispersion 
stratified  model  are  quite  different  from  the  barotropic 
not  clear  at  present  how  the  barotropic  curves  make  the 
the  stratified  curves.  This  question  will  be  addressed 


in  future  work. 
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The  two-layer  model  of  Bane  (1980)  has  shown  that  the  location 
of  the  density  front  in  relation  to  the  topography  is  very  important 
in  determining  the  dispersion  properties  of  the  various  waves. 

It  is  expected  that  the  same  is  true  of  the  continuously  stratified 
case;  however,  this  has  not  yet  been  investigated. 

The  horizontal  resolution  of  the  numerical  method  employed  in  this 
model  was  not  quite  as  good  as  desired.  Much  of  the  complex  structure 
in  figures  5.2  through  5.7  resulted  from  contouring  only  a  few  mesh 
points,  and  should  therefore  be  viewed  as  merely  suggestive  of  the 
true  wave  structure.  A  new  numerical  method  has  been  recently  developed 
at  NCAR  that  should  improve  the  resolution  of  this  model  by  an  order 
of  magnitude  in  future  computations. 

This  study  has  only  considered  small  amplitude,  stable  perturba¬ 
tions  to  the  basic  state.  The  problems  of  baroclinic  instability  and 
critical  levels  in  the  mean  flow  have  been  excluded  in  order  to  simplify 
the  analysis.  Inspection  of  the  quasigeostrophic  potential  vorticity 
field  (figure  5.8)  reveals  that  the  potential  vorticity  gradient  changes 
sign  seaward  of  the  frontal  cone.  This  is  a  necessary,  but  not  a  suf¬ 
ficient,  condition  for  baroclinic  instability  (cf  Pedlosky,  1979,  p.  440); 
therefore,  unstable  waves  are  possible,  but  a  more  detailed  analysis  is 
required  to  determine  the  exact  form  of  the  perturbations  for  which  the 
flow  is  unstable. 

Recent  measurements  in  the  Gulf  Stream  south  of  Cape  Hatteras  in¬ 
dicate  that  a  dominant  mode  of  oscillation  of  the  Stream  corresponds 
to  a  downstream  propagating  disturbance  with  a  period  of  eight  days 
(Bane  and  Rrooks,  1979b;  Brooks  and  Banc,  1980).  The  present  model. 
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however,  only  considers  onstream  propagating  disturbances,  and  thus 
avoids  the  possibility  of  critical  levels  in  the  flow.  A  logical 
extension  of  this  model  would  be  to  consider  downstream  propagating 
waves  and  unstable  waves  hv  allowing  the  wave  frequency  to  assume 
negative  and  complex  values. 

The  verical  wall  imposed  at  the  coast  is  not  a  very  realistic 
boundary  condition,  '['he  Kelvin  wave  mode  (mode  01  is  uurclv  a  conse- 
uucncc  of  this  boundary  condition.  The  vertical  wail  requires  a  node 
in  the  cross-shelf  velocity  at  the  coast.  A  more  realistic  sloping 
beach  boundary  condition  would  only  require  that  the  solution  remain 
bounded  at  the  coast;  however,  this  would  be  much  more  difficult  to 
treat  numerically.  Bane  (1980)  compared  the  sloping  beach  and  vertical 
wall  boundary  conditions  in  the  two- layer  case  and  found  that  the  ver¬ 
tical  wall  did  not  affect  the  solutions  appreciably  so  long  as  the  wall 
was  of  small  height  compared  to  the  total  depth  and  was  far  removed 


from  the  frontal  zone. 


7.  Conclusions 


It  has  been  shown  that  the  presence  of  a  continuously  stratified 
boundary  current  such  as  the  Gulf  Stream  in  a  coastal -trapped  wave 
model  permits  the  existence  of  a  new  class  of  fronta 1 -trapped  waves  in 
addition  to  the  topographic  Rossby  wave  (continental  shelf  wave)  modes. 
The  fronta 1 -t rapped  wave  is  a  Rossby-like  wave  in  that  it  is  supported 
by  the  ambient  quasi -gcostrophic  potential  vorticity  gradient  in  the 
cyclonic  side  of  the  mean  current.  This  potential  vorticity  gradient 
is  due  to  both  the  shear  of  the  mean  current  and  the  associated  sloping 
density  surfaces.  Banc  (1980)  also  has  found  frontal -trapped  waves  in 
a  two-lavcr  Gulf  Stream  model.  In  the  present  model,  the  mean  current 
flows  along  the  continental  slope,  with  its  inshore  surface  front  loca¬ 
ted  approximately  over  the  shelf  break.  The  wide  continental  shelf  sim¬ 
ulates  the  topography  off  Onslow  Bay,  North  Carolina.  The  stratifica¬ 
tion  over  the  shelf  is  level,  so  that  only  topographic  Rossby  waves 
may  exist  there.  front  a  1 -t rapped  waves  may  exist  in  the  frontal  cone 
of  the  mean  current.  These  two  regions  are  coupled,  and  the  resulting 
wave  motion  may  be  cither  a  shelf  wave,  a  frontal -trapped  wave,  or  a 
hybrid  between  the  two.  Following  a  particular  dispersion  curve,  the 
modal  structure  changes  from  one  wave  form  to  another,  indicating  a  mode¬ 
coupling  resonance  similar  to  that  found  by  Bane  (1980).  The  coupling 
between  the  two  wave  guides  (the  topographic  potential  vorticity  gradient 
and  the  quasi -gcostrophic  potential  vorticity  gradient)  is  much  more 
subtli  in  the  continuously  stratified  case  than  in  the  two  layer  case. 

It  is  inappropriate  to  speak  of  a  particular  modal  structure  in  the 
continuous  case,  since  a  wave  mode  exhibits  different  structures  in 
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different  regions  of  the  dispersion  plane. 

The  importance  of  these  stable  Rossby-like  waves  to  the  mesoscalc 
variability  of  the  real  Gulf  Stream  is  not  known.  The  free  modes  found 
in  this  model  reveal  a  possible  mechanism  for  sub-inertial  frequency 
disturbances  propagating  along  the  inshore  side  of  the  Stream.  The 
complexity  of  the  solutions  in  this  highly  simplified  model  of  the 
Gulf  Stream  indicate  that  the  response  of  the  real  Stream  depends  upon 
a  number  of  interacting  factors. 
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Appendix 

A1 .  Separation  of  time  dependent  flow  from  mean  flow 
The  full  equations  of  motion  are 
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where 


a  ;>  a  3 

n  =  -  +  u  — .  +  v  — -  +  w  _ 

At  Dx  3y  3z 

q  -  ui  +  vj  +  wl; 

f  is  the  vertical  component  of  the  earth's  rotation,  which  will 
be  assumed  constant, 

e  is  the  horizontal  component  of  the  earth's  rotation,  which 
will  be  assumed  negligible. 


and  v  is  the  molecular  viscocity.  Henceforth,  the  fluid  will  be 
assumed  inviscid,  i.c.  v  =  0.  liquation  (A1.4)  can  be  separated  into 
two  equivalent  equations  because  the  fluid  is  incompressible: 


Do  =  0 


(A1 .5) 
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The  dependent  variables  consist  of  a  mean  or  time- independent  component 
and  a  fluctuating  or  time-dependent  part.  The  mean  velocity  is  in  the 
y-direction  only,  so  that  the  mean  velocity  components  in  the  x  and  z 
directions  are  zero: 


u  =  u\  v  =  v  +  v ' ,  w  =  w' 

P  *  PD  ♦  P  ♦  P*  .  P  =  P0  +  P  +  P' 

where  p  =  PQgz  is  the  hydrostatic  pressure. 

The  equations  of  motion  are  averaged  over  a  period  of  time  much 
greater  than  a  typical  time  scale  for  the  time- dependent  motions,  i.e. 


r 


r  any  dependent  variable. 


A  Boussinesq  approximation  is  adopted,  that  is,  (p,p';)<<D  so  that 

o 

p  is  replaced  by  pQ  everywhere  except  in  the  buoyancy  equation  (A1 .5) . 
The  mean  flow  equations  are  then 
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lu'u'  +  V'u '  +w'u'  I 

*-  x  y  z  J 

(A1.7) 

1  u '  v  +  v'  v'  +  w' V'  I 
*-  x  y  z  J 

(A1.8) 

g  +  1  u '  w'  +  v'  w'  +  w'  w'  1 

Lx  y  z  J 

(A1.9) 

The  terms  in  brackets  are  the  x,  y  and  z  components  of  the  Reynolds 
stress  and  will  be  denoted  as  R^X\  R^  and  R^ 


,  respectively.  The 
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buoyancy  is  defined  as 


h  =  -  SR =  h  +  h  +  b ' 


Po  ° 


and  the  non-hydrostatic  pressure  per  unit  density  as 


I’-i'o 

v  =  — - =  n  +  v 1 

P„ 


liquations  (Al .7)  through  (A1.9)  now  become 


fv  =  it  +  \V 
x 


(Al .10) 


0  =  n  +  R 

y 


(Al .11) 


0  =  7T  -  b  +  Rl 


(Al .12) 


Equations  ( A 1 . 10)  through  fAl.12)  are  added  to  (Al.l)  through  (A1.3) 
to  give  the  equations  of  motion  for  the  time  dependent  flow: 


u1  +  vu'  +  u'u'  +  v ' u '  +  w'u’  -  fv'  = 

t  y  x  y  z 


-7i'  + 


(Al .13) 


v'  +w'  +  u'v'  +v'v'  +w'v'  +  fu'  +vu'  +vw' 
t  y  x  y  z  x  z 


-■n'  +  Rl 

y 


(Al . 14) 


w'  <■  vw'  +■  u'w'  +  v'w'  +  w'w' 
t  v  x  y  z 


=  —  tt  '  +  b'  + 


(Al .15) 


u '  +  v '  +  w '  =0 
x  y  z 


(Al .16) 


2  2 

b'  +  vb'  +  M  u '  +  N  w1  +  u’b'  +  v'b'  +  w'b' 
t  y  x  y  z 


0 


(Al .17) 
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Equations  (A1.13)  through  (A1.17)  are  the  full  non-linear  equations 
for  the  wave  motion.  In  this  analysis,  the  non-linear  terms  and  the 
Reynolds  stresses  are  assumed  to  be  small  compared  to  the  other  terms 
in  the  equations.  This  can  he  verified  through  a  simple  scale  analysis 
(cf.  Brooks,  197S,  pp.  212-218).  A  hydrostatic  approximation  is  adopted, 
since  the  vertical  velocity  w  is  small,  and  its  derivatives  for  long 
waves  are  smaller  still  when  compared  to  the  pressure  and  buoyancy  terms. 
The  linearized,  hydrostatic  equations  for  the  wave  motion,  with  the 
primes  dropped,  arc 


II  +  VU  -  l'v  =  -TT 
t  y  x 


v  +  vv  +  (f+v  )  u  +  v  w  =  -n 
t  y  x  z  y 


0  =  -7T  +  b 

z 


u  +  v  +  w  =0 
x  y  z 


+  vb  +  M“u  +  N^w  =  0 


(A1.18) 

(A1.19) 

(A1 .20) 

CA1.21) 


( A 1 . 22) 


A2.  Formulation  of  the  governing  equation. 

The  equations  of  motion  for  the  mean  flow  can  be  combined  by  dif¬ 
ferentiating  (A1.10)  with  respect  to  z  and  (A1.12)  with  respect  to  x, 
and  subtracting  to  give  the  thermal  w'.id  relation 


f v  =  M'  , 


(A2  1) 


which  will  be  applied  numerous  times  in  this  section. 

Alongshore  propagating  disturbances  are  sought;  therefore  the 

i(o0t+*y) 

perturbation  variables  are  taken  to  have  the  form  e  as  their 

(y.t)  dependence.  The  equations  of  motion  (A1.18)  through  (A1.22) 
now  become 


iou  -  fv  =  -7i 


(A2 . 2) 


i  o  v  +  ( f+  v  )  u  +  v  w  =  -  iZv 
X  z 


(A2.3) 


0  =  -  n  +  h 
z 


u  +  i£v  +■  w  =0 
x  z 


C A2 .4) 
(A2.5) 


2  2 

iob  +  Nfu  +  N~w  =  0 


(A2.6) 


where  o  =  is  the  intrinsic  frequency. 


F.liminating  the  buoyancy  b  between  (A2.4)  and  (A2.t>)  gives 


Mu  *  N  w  =  -ion 

z 


(A2.7) 
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Using  this  result  to  eliminate  w  from  (A2.3)  gives 


oN^v- i  [N~  ( f+v  )  -  M^v  ]u  =  -£N2it  +  ov  tt 
1  X  z  z  z 


(A2.8) 


Using  (A2.2)  to  eliminate  v  in  (A2.8)  and  applying  the  thermal  wind 
relation  (A2.1)  gives 


aN2  ,  M2  if  . 

u  =  — —  (n  -  -Tr  7T  +  —  n) 
lA  X  n2  z  o 


(A2.9) 


where  A  =  (a2-o^) -M4  and  a2  =  f(f+vx) .  Solving  for  v  by  elimina¬ 

ting  u  between  (A2.2)  and  (A2.8)  gives 


v  =  — [(ATo2  -  M4)tt  -  o2M2tt  +  of£N2ir] 

T  A  i  X  Z 


(A2.10) 


The  vertical  velocity  w  is  found  by  using  equation  (A 2.7)  and  (A2.9) 


io,..2  ,2  2.  M^f  , 

w  =  -£[m  *x  -  (of  -  o  )v~ r 


(A2.ll) 


Substituting  (A2.9),  (A2.10)  and  (A2.ll)  into  the  continuity  equation 
(A2.5)  gives  the  governing  equation 

f-oN2  ,  M2  If 

l  A  x  n2  z  a  > }  x 


£  (N2o^  -  M4)ttx  -  o2M2nz  +  of£N2Ttj 

[i  <"2\  -  <«f  ■  o 


(A2.12) 
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After  carrying  through  the  differentiations,  rearranging, 
the  identities 


N2  =  M2  ,  (  a2  )  =  M2 

X  z  \  f  Iz  X 

„  l  .2  2  2alM2 

z  f  ’  z  f 

1,2  c2,  -  „ 

°x  =  f  (0f  -  f  >  =  Vx£ 


all  of  which  follow  from  the  thermal  wind  relation  (A2.1) 
equation  becomes 


The  discriminant  of  (A2.13)  is 

A 


r 


N 


,  2  2. 

(of-o  ) 


N 


and  using 


,  the  governing 


(A2.13) 


(A2.14) 


which  is  negative  for  the  density  profile  considered  in  this  study  as 
long  as  a q<  f ;  therefore  the  governing  equation  is  elliptic. 


A3.  Finite  difference  form  of  the  forced  equation 


The  governing  equation  (2.25)  has  the  general  form 

An  +  Bn  +  cn  +  on  +  En  +  fit  =  o  (A3.i) 

xx  xr,  r,r.  x  z, 

where  A,  B,  C,  D,  H,  and  F  are  functions  of  (x.r  ;  o^,  ().  If  we 
consider  n  to  consist  of  a  perturbation  pressure  it,  plus  a  constant 
forcing  term  t,  then  the  forced  equation  is 

Air  +  Bn  +  Cn  +  Dtt  +  Ett  +  Fir  =  -Ft  (A3. 2) 

xx  xr,  r.r.  x  r. 

This  forcing  term  is  equivalent  to  forcing  the  system  by 

T  ei(o0t+£y) 

and  is  used  only  to  search  for  the  resonance  response  of  the  system. 

The  usual  second-order  accurate,  centered  difference  analogues 
of  the  derivative  terms  are  used  for  the  finite  difference  form  of  (A3. 2) 


A.  .  (it  .  .  .  -  2tt  .  .  +  m .  .  . ) 

ij  i  +  l  ,J  t .1  i-i.J 

i 

A  x“ 


B.  . 


l*  .  . 

1 J 

(n .  ,  .  , 

i+l, 3+1 

-  IT.  . 

1-1  ,J+1 

-  IT 

4AxA  C 

C.  . 

__LL 

(it.  .  - 

i,J  +  l 

2lT  .  .  +  TT  . 

j-i} 

A  2r, 

IJ  1, 

D.  . 

F..  . 

_LL 

(n  .  .  - 

1  +  1  ,J 

IT  .  .  .  )  + 

1-1,  J 

2Ax 

2A£ 

i-l} 


+  F.  .  IT .  .  =  C. .  . 
U  1J  iJ 


(A3. 3) 
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where  A .  . 

ij 


=  A(x . 
1 


V°o’f) 


etc . ; 


G.  .  =  -F. .t; 
13  13 


Ax  =  ,  Ac  =  j~  ■  an^  1  and  ^  are  t'le  number  of 

grid  points  in  the  x  and  £  directions,  respectively. 


Some  rearrangement  leads  to 


-B.  . 

U  n 

ihL 

+  Lij) 

n.  .  . 

1,3  +  1 

B.  . 

+  y  r 

4  Ax  A:  'i 

-l.j  +  1  + 

2A  J 

4AxA  C  i  +  l ,3  +  1 

/A 

+(-H  - 

'Ax' 

J±)  I! 
2Ax  /  1 

-i.j  * 

[f,  - : 

t A.  . 

VAx“ 

+ 

a^2ZI  y 

*( 

2  + 

'  Ax 

D.  A 
J^L) 

2  Ax/ 

ni+l  ,3 

B.  . 

♦  _LL_ 

*“T 

ACJX 

F..  . 

_y 

)  n.  . 
J  1.) 

+  f  By  )  n 

4  Ax  A  C 

i  -  1  ,  j  -  1 

l  2  A  j 

-1  V4AXAC/  i+l, 

(A3. 4) 

This  is  a  nine-point  stencil  over  the  interior  ^"id  points,  that  is, 
for  i  =  2  to  1-1  and  3  =  2  to  J-l  (See  Figure  A3.1). 

This  equation  is  solved  using  the  FORTRAN  subroutine  EVP9G,  gener¬ 
ously  provided  by  Pat  Roache.  The  subroutine  utilizes  the  Error  Vector 
Propagation  (EVP)  algorithm  described  in  Roache  (1978).  Given  the  co¬ 
efficients  for  the  finite  difference  equation  (A3. 4)  and  the  coefficients 
for  the  boundary  conditions,  EVP9G  returns  the  solution  for  7t. 
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The  boundary  conditions  (2.26)  through  (2.28)  are  of  the  form 


A  .  M 
nb  x 


B  .  IT 
nb  r. 


+  C  7T 
nb 


=  0 


(A3. 5) 


I- or  the  bottom  and  top  boundary  conditions  (subscript  nb  =  1,2  respec¬ 
tively),  I  VP9G  uses  the  second  order  accurate  finite  difference  approx¬ 
imation  for  the  x  derivative  but  allows  only  a  first-order  approximation 
to  the  normal  gradient.  A  second-order  form  is  possible  but  would  re¬ 
quire  a  bi -tridiagonal  matrix  solution  to  start  the  march. 

At  the  inshore  boundary  (nb  =  3)  both  the  depth  and  the  stratifi¬ 
cation  arc  very  small;  therefore,  to  greatly  simplify  the  computations, 
it  is  taken  to  be  zero.  The  boundary  condition  then  becomes 


Again,  the  finite  difference  approximation  to  the  normal  gradient  is 
only  accurate  to  first  order. 

The  offshore  boundary  condition  (nb  =4)  is  a  simple  Dirichlet  con¬ 
dition  and  its  values  are  stored  in  .  when  FVP9G  is  called. 

I ,} 


A 4 .  Subroutine  LVP9C  (I NIT,  F,  Cl,  C 2,  C3,  C4 ,  C5 ,  C6,  C7 ,  C8 ,  C9,  CIO, 


IL,  JL,  I I.D,  NCRS ,  NBJ1 ,  NBJL ,  CTA,  CTB,  CTC ,  CBA, 
CRB,  CBC,  NBI1,  CLEF) 

This  subroutine  solves  the  equation 

Cl ( 1 ,J)*F(I-1 ,J+1) +C2( I , J) *F(I ,  J  +  l ) +C3 ( I , J) *F ( I  +  1  ,J  +  1) 

+C4  ( I , J) *F(  I- 1  , J) +C5 ( I , J) *F (I , J) +C6( I , J) *F( 1+1  ,J) 

♦C7(I,J)*F(I-1 ,J-1)+C8(I,J)*F(I,J-1)  +  C9(I,J)*F(I+1,J-1)  =  C10(1,J) 

by  marching,  using  the  LVP  algorithm  (see  Roache,  1978) 

Calling  sequence: 

CALL  l.VP9G(  I N  IT ,  F,  Cl,  C2,  C3,  C4 ,  C5 ,  C6 ,  C7 ,  C8 ,  C9,  CIO,  IL,  .JL 
I LD,  NCRS,  NBJ 1 ,  NBJL,  CTA,  CTB,  CTC,  CBA,  CBB ,  CBC, 
NB11,  CLEF) 

Pa  raniet  ers : 

INIT:  Flag  to  determine  whether  the  call  is  a  repeat  solution  (IN1T=0) 

or  an  initialization  call  ( I NIT= 1 ) . 

F:  Matrix  of  dimension  (ILD,  JL) .  On  output,  F  contains  the  solution 
to  the  finite  difference  equation.  If  Dirichlet  conditions  are 
used  at  any  boundary,  the  boundary  values  are  set  by  whatever  values 
are  stored  in  the  appropriate  position  in  F  when  the  subroutine  is 
called,  e.g.  for  the  offshore  boundary,  the  values  stored  in 
I  (1I..J)  .  J  -  1,  JL  arc  used  as  the  boundary  values. 

Cl,  (12 ,  .  .  .  CIO:  Matrices  of  dimension  (ILD,  JL)  .  On  input  they 


contain  the  coefficients  of  the  finite  difference  equation  as 
given  above.  The  values  are  not  changed. 
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II.:  Number  of  grid  points  in  the  x-direction.  Must  be  less  than  or 

equal  to  ll.P. 

Number  of  grid  points  in  the  y-direction.  Must  be  less  than  or 
equal  to  ILD. 

ll.P:  first  dimension  of  F  in  calling  program.  In  this  version  of 
FVP'.Ki ,  ILD  must  be  less  than  or  equal  to  31,  but  this  can  be 
changed  by  changing  the  appropriate  dimension  statements. 

NCRS:  Number  of  cross  derivatives.  If  Cl  =  C3  =  C7  =  C9  =  0,  setting 
NCRS  =  0  saves  time  in  the  computation.  Otherwise,  NCRS  =  1, 
and  the  subroutine  GTRI ,  a  tridiagonal  solver,  must  be  made 
available  to  FVP9G.  GTRI  solves  for  the  (J+l)th  row  in  the  march. 

N'R.J  1 :  Flag  to  determine  the  form  of  the  boundary  condition  at  the  bottom 
(./  =  I)  boundary.  NBJ1  =  1  gives  the  general  oblique  derivative 
form  of  the  boundary  condition  CBA( I ) * ( F ( 1 + 1 , 1) - F ( I  - 1 , 1) ) +CBB ( I ) * 
(F( I , 2 ) - F ( I , 1 ) )  +  CBC( I ) *F ( I , 1)  =  0. 

NBJl=d  gives  simple  Pirichlet  conditions.  In  this  case,  the 
boundary  values  are  stored  in  F(I,1),  I  =  1,  IL  on  input  and  CBA, 
CRB  and  CRC  are  ignored. 

NB.JL:  Flag  to  determine  the  form  of  the  boundary  condition  at  the  top 
.  Analogous  to  NBJ  1  . 

CIA,  (IB,  CTC :  If  NB.JL  =  1,  arrays  of  length  ILD  containing  the  coef¬ 
ficients  for  the  oblique  derivative  boundary  condition  at  the 
top  (,I=JL)  : 

CTA  ( 1 )  *  (  F  ( I  + 1  ,.JL)  -F(I- 1  ,.JL)+CTB(I)  *(F(I,JL)-F(1  ,JL-1)) 


+CTC ( 1)*F(I ,JL)  =  0 
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CRA ,CBB ,CBC :  If  NBJ1  =  1,  arrays  of  length  ILD  containing  the  coeffi¬ 
cients  for  the  oblique  derivative  boundary  condition  at  the 
bottom  (.1=1).  See  NBJ 1  above. 

NBI1:  Flag  to  determine  form  of  boundary  condition  at  the  left  boundary 
(1  =  1).  For  a  Dirichlet  condition,  NB 1 1  =  1 ,  and  the  boundary 
values  arc  stored  in  F(1,J),  J=1 ,  JL.  For  a  mixed  condition 
(sometimes  called  a  Robin's  condition),  i.e. 

F  /F=C 
x 

NB  11=3.  This  version  of  F.VP9G  does  not  allow  oblique  deriva¬ 
tive  conditions  at  the  left  boundary. 

CLFF:  If  NBI1=5,  it  is  the  constant  in  the  finite  difference  form  of 

the  mixed  boundary  condition  at  the  left  boundary: 

( F(2 , J) - F( 1 , J) ) / F ( 1 , J) =CLEF .  Otherwise  it  is  ignored. 

Precision:  All  arrays  in  the  calling  sequence  are  double  precision. 

Required  supporting  subroutines:  F.VP9G  requires  the  double  precision 
version  of  UNPACK  matrix  algebra  program  along  with  the  BLAS 
(basic  linear  algebra  system)  which  LINPACK  uses.  If  NCRS=1 , 
the  tridiagonal  solver  GTR1  is  also  required.  GTRI  is  provided 
with  this  version  of  EVP9G. 

Output:  The  solution  matrix  is  returned  in  F.  If  the  system  is  un¬ 
stable,  an  error  message  is  printed.  On  repeat  calls  (INIT=0) 
the  maximum  error  and  the  final  errors  are  printed. 

Programming  notes: 

1)  If  the  oblique  derivative  form  of  the  boundary  condition  is  used 

and  either  CTB  or  CBB  are  zero,  that  is,  the  normal  gradient  is  not 


involved,  a  zero-divide  results. 
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2)  If  CLF.F=1.0,  a  zero  divide  will  result. 

3)  The  only  boundary  condition  allowed  at  the  right  boundary  is  a 
simple  Oirichlet  condition. 

4)  When  there  is  no  cross-derivative  term,  it  is  not  necessary  to 
set  NCRS=0,  but  it  saves  some  time. 
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Fig.  2.1: 
Fig.  3.1: 

Fig.  3.2: 


Fig.  3.3: 

Fig.  3.4: 

Fig.  4.1: 

Fig.  4.2: 

Fig.  4.3: 


Schematic  of  idealized  model  geometry. 

Map  of  the  southeastern  United  States  coast,  showing  depth 
contours  and  average  position  of  the  Gulf  Stream. 

Average  Gulf  Stream  velocity  section  taken  off  Cape  Fear 
in  summer,  1968.  Contours  are  isotachs  of  downstream 
velocity  in  cm/ sec .  Taken  from  Richardson,  Scmitz  and 
Niiler  (1969)* 

Model  mean  velocity  field  and  bottom  topography.  Contours 
are  isotachs  of  downstream  velocity  in  cm/sec.  Vertical 
and  offshore  coordinates  are  dimensionless. 

Model  mean  temperature  field.  Contours  are  in  degrees  Cel¬ 
sius  . 

Solution  grid.  Crosses  mark  the  grid  points  used  in  the 
finite  difference  approximation. 

Integrated  kinetic  energy  versus  dimensionless  frequency 
for  Z  =  2.0.  The  eigenvalues  of  oq  are  found  by  locating 
the  local  maxima  of  the  kinetic  energy. 

Eigenvalue  curves  for  the  test  case.  The  curves  labeled 
I  and  II  correspond  quite  well  with  the  Mode  I  and  Mode  II 
topographic  Rossby  wave  dispersion  curves  of  Wang  and  Mooers 
(1976),  which  are  shown  in  the  inset  as  solid  lines.  The 
present  model  was  able  to  resolve  several  other  modes,  in 
addition  to  those  identified  by  Wang  and  Mooers  (1976) . 
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Fly.  5.1:  Dispersion  curves  obtained  from  the  model  geometry  shown 
in  Fig.  3.3  and  3.4. 

Fig.  5.2:  Structure  of  the  mode  1  wave  for  small  wavenumber  [l  -  1.25). 

(a)  Perturbation  pressure  ti  .  Note  the  barotropic  portion 
over  the  shelf  transitioning  to  a  strongly  baroclinic  struc¬ 
ture  in  the  inshore  side  of  the  mean  current .  Contours  are 
in  dimensionless  units.  (b)  Alongshore  velocity  v,  and 
(c)  cross-shelf  velocity  u.  Note  again  the  strong  baroclin- 
icity  in  the  cyclonic  side  of  the  current  transitioning  to  a 
barotropic  structure  over  the  shelf.  The  velocity  contours 
roughly  follow  constant  temperature  (density)  surfaces  (see 
Fig.  3.4). 

Fig.  5.3:  Structure  of  the  mode  1  wave  for  intermediate  wavenumber 
(£  =  2.0).  Shown  here  is  the  perturbation  pressure  it  in 
dimensionless  units.  The  barotropic  component  over  the 
shelf  is  absent  at  this  wavenumber.  The  mode  1  wave  is  a 
pure  frontal -trapped  wave  at  this  wavenumber,  with  its  struc¬ 
ture  confined  to  the  surface  layers  on  the  cyclonic  side  of 
the  mean  current . 

Fig.  5.4:  Perturbation  pressure  n  structure  for  the  mode  1  wave  for 
large  wavenumber  (£  =  4.0).  The  frontal -trapped  structure 
has  disappeared,  and  this  mode  has  become  a  barotropic  continental 
shelf  wave,  or  topographic  Rossby  wave. 

Fig.  5.5:  Structure  of  the  mode  2  wave  for  intermediate  wavenumber 
{l  =  2.0).  (a)  Perturbation  pressure  it.  Note  the  slight 

bottom  trapping  due  to  the  density  stratification,  (b)  Along¬ 
shore  velocity  v.  (c)  Cross-shelf  velocity  u. 
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Pig.  5.6: 

Fig.  5.7: 
Fig.  5.8: 
Fig.  5.9: 

Fig.  6.1: 

Fig.  A3. 1 : 


Perturbation  pressure  n  for  the  mode  2  wave  at  high 
wavenumber  (?=  4.0). 

Perturbation  pressure  it  structure  for  the  mode  3  wave. 
Perturbation  pressure  tt  structure  for  the  mode  4  wave. 
Quasigeostrophic  potential  vorticity  for  the  undisturbed 
mean  current. 

Barotropic  shelf  wave  dispersion  curves  for  the  model  bottom 
topography  with  no  mean  current  or  stratification. 

(a)  Nine  point  finite  difference  stencil  used  in  solving 
the  governing  equation,  (b)  Schematic  of  solution  mesh. 
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